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Measurements of QCD observables now have reached such high precision that power 
corrections to the structure functions can often be extracted with a reasonable accuracy 
from the existing data. The situation on the theoretical side is much less clear. In 
the best understood case of deep inelastic scattering, the relevant contributions can be 
attributed in the framework of operator product expansion (OPE) to matrix elements 
of higher twist operators [|l|, but their determination in QCD is ambiguous due to 
occurrence of power divergences [0]. From the phenomenological point of view, however, 
attempts to compute these matrix elements using e.g., QCD sum rules or the bag model 
have provided results which seem to have at least the right order of mag nitude i, i |, i 
as compared with available experimental estimates. 

Recently there has been some interest in another phenomenological approach to 
power-suppressed corrections in QCD |0, based on the fact that the only possibility to 
interpret the higher-order radiative corrections in a consistent manner (i.e. as asymptotic 
series) requires the existence of power suppressed terms The IR-renormalon contri- 
butions occur because certain classes of higher order radiative corrections to twist-2 are 
sensitive to large distances, contrary to the spirit of OPE. Although these divergences 
have to cancel with UV power divergencies in matrix elements of twist-4 operators i.e., 
they are totally spurious, there are arguments which suggest that they reflect the cor- 
rect shape, if not the magnitude, of the power-suppressed terms |T0|, |ll[]. In this sense 
we define the prediction of the renormalon model of higher-twist corrections to a DIS 
structure function as the power-suppressed uncertainty which occurs in the perturbative 
expansion of the Wilson coefficient to the twist-2 contribution. 

Commonly the divergent series of radiative corrections is regarded as an asymptotic 
series and defined by its Borel integral. The actual calculations are done in a large- iVp 
limit which allows to resum the fermion bubble-chain to all orders yielding the coefficient 
of the a^Np~'^ - term exactly. Subsequently, it is converted into the exact coefficient of 
the ag(3Q~^ - term by the substitution Np — > Np — 33/2 = — 67r/3o, known as the 'Naive 
Non-Abelianization' (NNA) (see the last two references of ||^). The asymptotic character 
of the resulting perturbative series leads to resummation ambiguities - a singularity in 
the Borel integral destroys the unambiguous reconstruction of the series and shows up 
as a factorial divergence of the coefficients of the perturbative expansion. The general 
uncertainty in the perturbative prediction can be estimated to be of the order of the 
minimal term in the expansion, or by taking the imaginary part (divided by vr) of 
the Borel integral. Both procedures lead to resummation ambiguities of the form C x 
{A? /Q'^Y , with r = 1 for the leading IR renormalon, and a numerical coefficient C which 
either can be taken as it comes out from the NNA calculation, see below, or can be fitted 
as a free parameter. 

The unpolarized hadronic scattering tensor for leptons scattering off nucleons can be 
divided into three structure functions ||12| 
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Here x = Q'^/{2p ■ q) and Q"^ = —q^. The structure functions Fi,i = L,2,3 can be 
generally decomposed in the following manner 

F,{x, Q') = Ft\x, Q') + Q') + ^,K{x, Q') + , (2) 

where describes the leading twist-2 contribution, /i™^ describes the target mass 
corrections which are directly related to twist-2 matrix elements [|13|. hi contains the 



genuine twist-4 contribution which is in principle sensitive to multiparton correlations 
within the hadron, and which we want to estimate using the renormalon model. Note 
that hi has dimension 2 (if radiative corrections are neglected) and hence it has to be 
proportional to a certain mass scale /i^. In the analysis of the experimental data used in 
this paper the coefficients hi were extracted by accounting for target mass corrections up 
to order ^(M^/Q^), except for the case of /il where target mass corrections have been 
accounted for to the order 0{Mjf/Q'^) only. In the analysis of the twist-4 contributions 



T^, |T5|, [T^ it is common to neglect the Q dependence of hi, which is due to radiative 
corrections, and to sum target mass corrections and the twist-2 contributions into a 
leading-twist (LT) structure function thus arriving at the notation 

F,{x,Q') = Fr\x,Q') + ^^h,{x) = Fr\x,Q') U + ^) • (3) 



Now we shall shortly summarize the main features of the model. In the renormalon 
approach the expression for the higher twist correction hi{x) to a structure function Fi 



has the form of a Mellin convolution 17 



hi = Pi0 Ft^ (4) 

of a coefficient Pi{z), and the twist-2 part Fl~'^{x) of the structure function Fi. The mass 
scale ix^ which enters Pi{z) equals to = SnCpA'^e''-^'' / Pq where Cp = 4/3, As is the 
QCD scale parameter in a given renormalization scheme s, and Cg is determined by the 
finite part of the fermion loop, Cjj^ = —5/3. This identification, which follows from the 
NNA prescription, is RG invariant only in the strict large-A^j limit. Another difficulty 
stems from the observation that within the logic of the model one cannot distinguish 
between LO and NLO parametrizations of twist-2 structure functions Fl~'^{x). In the 
following we will consequently use the LO GRV parametrizations of parton distributions 
but we have checked that qualitatively the same conclusions follow from a calcu- 



lation which incorporates the NLO parametrizations. In the present analysis we always 
refer to a region of around 2 GeV^ i.e., below the charm threshold, and accordingly 
we adopt the A^f = 3 value of the A parameter from the LO GRV fit, A = 232 MeV. We 
mention also that due to the two possible contour deformations above or below the pole 
in the Borel integral, as discussed above, in principle the overall sign of the coefficients 
Pi{z) cannot be uniquely determined. 
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We shall state from the very beginning that because the IR-renormalon result is 
proportional to the twist-2 contribution, it cannot be expected to describe the com- 
plete twist-4 correction which contains a genuine multiple field correlation, and which 
depends therefore on the exact internal hadron wave function. This distinction has a 
real physical meaning as can be seen from the following argument. Let us assume that 
the same structure functions F{x, Q"^) are measured for different hadrons and that the 
lowest-order (in a^) leading twist and l/Q^ parts {F^~^{x)) can be separated 

experimentally. Then the renormalon contribution cancels for the difference of the ratio 
of moments 

which is, however, still sensitive to the genuine higher-twist corrections, which depend on 
the exact specific internal hadron wave function. We hope that this argument elucidates 
the fundamental limitations of the renormalon approach. Nevertheless, to achieve a 
better understanding of the renormalon contributions it is important to study as many 
different CcLSGS cLS possible. On the basis of these results we hope to develop a physical 
interpretation which explains why this phenomenological approach is successful in some 
cases, but fails in others [|T^ . 

So far the renormalon model was used to calculate the power corrections for the 



non-singlet part of the structure functions -^l ll^l ^^^1 gi |20|. The physically equivalent 
approach based on dispersion relations of ^ was used to calculate the higher twist 
contribution to the non-singlet part of the structure functions F2 and F3 p2 . 

In the case of the structure function F2{x) an experimental analysis of higher twist 



corrections exists for proton and deuteron targets |T^. Assuming to be low enough, 
i.e., below the charm threshold, one obtains 
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where -F2 i = 0, 3, 8 denote corresponding SU(3) combinations of parton densities, and 
G is the gluon twist-2 structure function of the nucleon. The corresponding expression 
for the Fy,{x) structure function reads 



(7) 



As already mentioned, so far only the coefficients Pns and P3 are known. For a first try, 
we set in (|]) the coefficient Pg to zero, and approximate Ps ~ Pns, which results in 
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Equations (|^), (J^), and define the renormalon model of higher twist corrections. 
In the comparison of the model predictions with the data we have used the NMC analysis 
of F2{x, Q"^) for C2,proton and C2,deuteron as Well as the more precise one for C2,proton - 
C2,neutron- For F3 wc take the new fit of [|16[. We have also confronted our results with 



the output of the computer code kindly provided to us by B. Webber 12^, but with the 



LO GRV [|T8| parametrization consistently used to represent twist-2 structure functions. 



Fig. |I| shows our results for Cp{x) and Cd{x) of F2, with twist-2 structure functions 



parametrized according to the leading order GRV-parametrization ||T8[, as well as the 
experimental results. Obviously the renormalon model underestimates the experimental 
data. On the other hand the proton and deuteron data differ only shghtly, indicating 
the dominance of the flavor singlet contribution. 

The x-dependence of the higher twist correction can be, however, reproduced very 
well if one adopts the value of /i^ a factor of 3 larger than what follows from the renor- 



malon model, /i^ as it was done in Ref. ||2^. Strictly speaking, the authors of 

used twist-2 MRS (A) ||23| parton distributions normalized at the scale = 10 GeV^ 
and hence they obtained /x^ ~ 2A^\. Note that as this scale dependence is governed by 
anomalous dimensions of twist-2 rather than twist-4 operators, it differs from what one 
would expect in QCD. 

It is interesting to speculate how taking into account the yet unknown coefficients Ps 
and Pq may influence the renormalon model predictions. First, note that the character- 
istic rise of the renormalon model prediction in Figure |^, which clearly follows the trend 
seen in the data, results from the term in Pns proportional to in the Mellin space 



20| , |2^ . If one assumes that Ps and Pg follow the same pattern, i.e. they do not contain 
terms proportional e.g. to iV^, then the gluonic contribution should be always negligible 
at large x. Indeed, the twist-2 gluon distribution should fall at least one power of 1 — x 



faster than the quark one |24|, and therefore it cannot be responsible for the discrepancy 
observed in Figure |I[ One could then try to constrain the singlet coefficient Ps taking 
into account the experimental fact that Cd{x) ~ Cp(x), and known relations between 
various parton distributions when x 1 ||25|, but it leads to the relation Ps = P^s? 
which has been discussed already above. Finally, we note that the deuteron and proton 
contributions to Eq. (|^) contain one unknown combination |Ps ® P2°'*^^ + Pq^G which 
therefore can be determined independently from Cp and Cd data: 

^Ps ® P2°'*"' + PG^G = Cd F^-'-' - -^Pns ® P|'*"' , (9) 



and 



^Ps ® P"''"' + Pg ® G = Cp Fi''-' - Ipns ® (P^'*"' + ^P^'*" 

9 DO 



We have explicitly checked that both equations and ( ]T0| ) are indeed consistent with 
each other, see Figure |^, and so one can take e.g. the arithmetic average between both 
estimates as the approximation for the unknown flavor-singlet contribution. Note that 
the curve is negative for x < 0.35 in accordance with the fact that the renormalon 
model prediction for Cp and Cd is higher than the data in this region, see Figure 
It turns out, however, that when one wants to use this approximation for comparison 
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with the difference C2,proton — C2,neutron-, which has been extracted in [jT4| with much 
smaller experimental errors, one encounters a problem. The results are shown in Figure 
^ together with the now fixed renormalon model prediction (including the fitted quark 
singlet and gluonic parts). The agreement with the data becomes better, but the model 
still underestimates the higher twist corrections for x between 0.2 and 0.4. 

In addition we have compared the renormalon model prediction for the difference 
of coefficients {h2p{x) — h2n{x))/x as given in Ref. |]T3|, see Figure |. As far as the 



x-dependence is concerned, the general trend seems to be reproduced by the model, but 
the scale has again to be taken a factor of 3 or larger than the We note that as 
{h2p{x) — h2n{x))/x is free from fiavor-singlet contributions, it is an ideal observable for 
models of higher-twist corrections, and that one should try to get better data (smaller 
error bars) for this quantity. 

The case of the purely non-singlet structure function F^{x) is illustrated on Figure 
^. Again, although the renormalon model prediction seems to have the correct shape, it 
falls short by a factor of 2 or 3 with respect to the magnitude of the higher-twist effects, 
as required by the data according to the LO analysis of \W\. Note that while the analysis 
of the magnitude of the coefficient xh^{x) in |16[ shows that it slightly decreases from 
LO to NNLO fits, its characteristic shape, indicating an increase from a negative to a 
positive value somewhere at 0.2 < a; < 0.65, remains stable. Finally, an analysis along 



the lines of Ref. shows that the renormalon model underestimates the magnitude 
of higher-twists corrections to by a factor of 2, see figure ^ Adjusting the scale /i^ 



to the higher twist corrections to F2{x,Q'^), as in p2|, results this time in a prediction 



which is larger than what follows from the present data. 

Summarizing, we conclude that the model in its present form can certainly be con- 
sidered as a useful phenomenological guide for estimates of higher-twist contributions 
to DIS, but one has still to be careful in claiming a phenomenological success of this 
approach. Our analysis shows that although the renormalon model does in fact a good 
job as far as the x-dependence of the higher-twist correction is concerned, the mass scale 
it predicts is by a factor 2-3 too low, depending on the observable. In other words, 
the NNA estimates of the leading powered-suppressed corrections to DIS turn out to be 
a factor 2 or 3 smaller then what seems to be required by the data. Other estimates. 



for example the results presented in |2^, should be seen from the same perspective. 
Another possibility would be to fit the mass scale independently to each observable, a 
point of view which has been recently advocated in |TI[]. One should also keep in mind 
that the analysis of ^2(0;) data still reveals some discrepancies and it is difficult to say 
now whether taking into account the singlet part, possibly with its specific mass scale, 
will improve the situation. We hope that these observations will help to understand 
successes and limitations of the renormalon phenomenology. 



Acknowledgment: The authors are grateful to A. Briill and B. Webber for many useful 
discussions. This work has been supported by BMBF, DFG (G. Hess Programm), DESY, 
and German-Polish exchange program X081.91. A.S. thanks the MPI fiir Kernphysik in 
Heidelberg for support. L. M. was supported in part by the KBN grant 2 P03B 065 10. 



6 



References 



[1] E.V. Shuryak and A.I. Vainshtein, Nucl. Phys. B199 451 (1982); 
R.L. Jaffe and M. Soldate, Phys. Rev. D26 49 (1982); 
R.K. Ellis, W. Furmanski, and R. Petronzio, Nucl. Phys. B207 1 (1982). 

[2] G. Martinelh and C.T. Sachrajda, Nucl. Phys. B478 660 (1996). 

[3] V.M. Braun and A.V. Kolesnichenko,Nucl. Phys. 283B, 723 (1987); 

I.I. Balitsky, V.M. Braun, and A.V. Kolesnichenko, Phys. Lett. B242 245 (1990); 
Phys. Lett. B318 648 (E) (1993). 

[4] E. Stein, P. Gornicki, L. Mankiewicz, A. Schafer and W. Greiner, Phys. Lett. B343 
369 (1995); 

E. Stein, P. Gornicki, L. Mankiewicz, and A. Schafer, Phys. Lett. B353 107 (1995). 

[5] V.M. Braun, QCD Renormalons and Higher Twist Effects, talk given at 30th Ren- 
contres de Moriond: QCD and High Energy Hadronic Interactions, Meribel les 
Allues, France, 19-25 March 1995, published in the Proceedings, Piep-ph/9505317l 
(1995). 

[6] X. Ji and P. Unrau, Phys. Lett. B333 228 (1994); 
M. Stratmann, Z. Phys. C60 763 (1993); 

X. Song and J.S. McCarthy, Phys. Rev. D49 3169 (1994), Erratum ibid. D50 4718 
(1994). 

[7] M. Neubert, Phys.Rev. D51 5924 (1995); 

M. Beneke, V.M. Braun, Nucl. Phys. B454 253 (1995); 
R. Akhoury, V.I. Zakharov, Phys. Lett. B357 646 (1995); 

C. N. Lovett-Turner, C.J. Maxwell, Nucl. Phys. B452 188 (1995) ; 

N. V. Krasnikov, A. A. Pivovarov, Mod. Phys. Lett. All 835 (1996); 
P. Nason, M.H. Seymour, Nucl. Phys. B454 291 (1995); 
P. Korchemsky, G. Sterman, Nucl. Phys. B437 415 (1995); 

M. Dasgupta, B.R. Webber, Power corrections and renormalons in e~^e~ fragmen- 
tation, CAVENDISH-HEP-96-9, |hep-ph/96085^ . 

[8] G. 't Hooft, in "The Whys of Subnuclear Physics," Erice 1977, ed. A.Zichichi 
(Plenum, New York 1977); 

A.H. Mueller, The QCD perturbation series., in "QCD - Twenty Years Later", 
edited by P.M. Zerwas and H.A. Kastrup, 162, (World Scientific 1992); 
V. I. Zakharov, Nucl. Phys. B385 452 (1992); 
A.H. Mueller, Phys. Lett. B308 355 (1993); 
M. Beneke, Nucl. Phys. B405 424 (1993); 

D. J. Broadhurst, Z. Phys. C58 339 (1993). 

[9] Yu. L. Dokshitser, G. Marchesini, B. R. Webber, Nucl. Phys. B469 93 (1996). 



7 



[10] Renormalons and Power Corrections in QCD, Proceedings of Nordita Worskhop, 
Copenhagen, Denmark, August 2-4, 1996, V. Braun and L. Magnea (Eds.), 
NORDITA-96-64-P, September 1996. 

[11] M. Beneke, V.M. Braun, L. Magnea, Phenomenology of Power Corrections in Frag- 
mentation Processes in e~^e~ Annihilation, |hep-ph/9701309| (1997). 

[12] A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta Phys. Rev. D18 3998 (1978). 

[13] W. R. Frazer, J. F. Gunion, Phys. Rev. Lett. B45 1138 (1980). 

[14] The New Muon Collaboration, Nucl. Phys. B371 3 (1992). 

[15] M. Virchaux and A. Milsztajn, Phys. Lett. B274 221 (1992). 

[16] A. V. Sidorov, Next to next-to-leading order QCD analysis of combined data for xF^ 
structure function and higher twist contribution. Presented at 13th International 
Seminar on High-energy Physics Problems (ISHEPP 96), Dubna, Russia, 2-7 Sep 
1996. |hep-ph/ 96093451 (1996) ; 
A. V. Sidorov, Phys. Lett. B389 379 (1996). 

[17] E. Stein, M. Meyer- Hermann, L. Mankiewicz, and A. Schafer, Phys. Lett. B376 
177 (1996). 

[18] M. Gliick, E. Reya, A. Vogt, Z. Phys. C67 433 (1995). 

[19] R. Akhoury and V.I. Zakharov, A Quick Look at Renormalons, piep-ph/9610492| 
(1996); 

R. Akhoury and V.I. Zakharov, Renormalon Variety in Deep Inelastic Scattering, 
|hep-ph/9701378| (1997). 

[20] M. Meyer-Hermann, M. Maul L. Mankiewicz, E. Stein, A. Schafer, 
Phys. Lett. B383 463 (1996), and Erratum to be pubhshed. 

[21] M. Beneke, V.M. Braun, Phys. Lett. B348, 513 (1995); 

P. Ball, M. Beneke, V. M. Braun, Nucl. Phys. B452 563 (1995). 

[22] M. Dasgupta, B. R. Webber, Phys. Lett. B382 273 (1996). 

[23] A. D. Martin, W. J. Stirling, R. G. Roberts, Phys. Rev. D50 6734 (1994). 

[24] S.J. Brodsky, M. Burkardt and I. Schmidt, Nucl. Phys. B441 197 (1995). 

[25] W. Melnitchouk, A.W. Thomas, Phys. Lett. B377 11 (1996). 

[26] L.S. Whitlow, S. Rock, A. Bodek, S. Dasu and E.M. Riodan, Phys. Lett. B250 193 
(1990). 

[27] The New Muon Collaboration, Phys. Lett. B364 107 (1995). 



8 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0. 

X 



Figure 1: Renormalon contribution to the coefficient Cp and Cd according to Eq. (^. 
The sohd hne shows the renormalon model predictions for Cp and C^, which are indis- 
tinguishable, obtained with the LO GRV parametrization |TB[. The dashed line shows 
the fit, as in Ref. |2^, which corresponds to the same x dependence, but the value of the 
scale parameter /x^ about 3 times larger than the one which follows from the renormalon 
model 



The filled and empty circles display the data for Cp and Cd according to Ref. 
respectively. 
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Figure 2: Determination of the unknown flavor- singlet component of the renormalon 
model from Cp (filled circles) and Cd (empty circles), see Eqs. (^) and (p!0|) respectively. 
The LO GRV parametrizations [|18] normalized at = 2 GeV^ have been used to 
generate the twist-2 parton distributions. The solid line is the arithmetic average. 
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Figure 3: The renormalon model prediction for Cp — Cn using the LO GRV parametriza- 



tion 



sohd hne, where C„ and Cd have been calculated according to Eq. 



Note 



that due to its definition, see Eq. Cp — Cn, is not a pure non-singlet quantity. The 
dashed line shows the prediction for Cp — Cn after the unknown fiavor-singlet contri- 
bution has been adjusted to reproduce the data for Cp and C^, see Eqs. (|^) and (|T0|). 
The dashed-dotted line shows the result of the same procedure, but with the scale /i^ 
adjusted, like in [^], to the description of the coefficients Cp and Cd, see Figure (|I]). 
The data points are from Ref. [|14| . 
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The renormalon model prediction for {h2p{x) — h2n{x))/x using the LO GRV 
], dashed hne. The dot-dashed hne shows the prediction with the 
scale /i^ adjusted to the description of the coefficients Cp and Cd, as in [Q. The data 
points are from Ref. Error bars for data points lying below x = 0.2 have been cut 



Figure 4: 
parametrization 



to fit into the figure. 
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Figure 5: Renormalon prediction for xh^lx) using the LO GRV [|I8| parametrization 
(solid line). The data points taken from Ref. |]l6l correspond to the LO analysis. The 
dashed line shows the prediction with the scale //^ adjusted to the description of the 
coefficients Cp and C^, as in [^. The dot-dashed line shows the original prediction of 
[|22| , obtained using the MRS A parametrization |]23[ normalized at = 10 GeV^. 
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Figure 6: Coefficient hip — h^n, obtained from tlie plienomeno logical parametrizations 
to R{x, Q"^) and F2 (solid line). The renormalon model prediction for h^p — h^n 
using the LO GRV parametrizations is depicted by the dashed line. The dot-dashed line 
shows the prediction with the scale [j? adjusted to the description of the coefficients Cp 



and Crf, as in 
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